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Abstract. - We identify a fundamental structure in the spectrum of microwave driven alkali 
Rydberg states, which highlights the remnants of the Coulomb symmetry in the presence of a 
non-hydrogenic core. Core-induced corrections with respect to the hydrogen spectrum can be 
accounted for by a perturbative approach. 



Introduction. - The excitation and subsequent ionization of Rydberg states of atomic 
hydrogen by microwave fields is one of the most prominent examples of the manifestation 
of classically nonlinear dynamics in a realistic physical system [Q. Given a driving field 
frequency comparable to the classical Kepler frequency of the unperturbed Rydberg electron, 
the electron's classical trajectory goes chaotic for sufficiently large driving field amplitudes, 
finally leading to its ionization on a finite time scale Correspondingly, large ionization 
rates are observed in experiments on real (i.e., quantum) Rydberg states of atomic hydrogen, 
in the appropriate parameter range |Q, ^ . 

As a matter of fact, already before the onset of classically chaotic motion, i.e. at not 
too large driving field amplitudes, individual quantum eigenstates of the atom in the field 
exhibit energies and ionization rates which are determined only by the orbital parameters 
of the classical trajectory they are associated with Those orbits which are the least 
stable under the external perturbation (i.e., which turn chaotic for the lowest values of the 
driving field amplitude, such as straight line orbits parallel to the field polarization axis for 
a linearly polarized drive) induce the largest ionization rates for their associated eigenstates. 
Consequently, in this near-integrable regime of classical dynamics, it is possible to classify 
the eigenstates of the atom in the field through quantum numbers associated with the orbital 
parameters of unperturbed Kepler ellipses, i.e. with the angular momentum and the Runge- 
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Lenz vector. An adiabatic invariant governs the slow evolution of these parameters under 
external driving Q. 

It should be noted, however, that a considerable part of experimental data has been 
accumulated in experiments on Rydberg states of alkali atoms rather than of atomic hydrogen 
[H 0' B 0' ^ pi'iori, a classical-quantum correspondence as briefly sketched above for 
atomic hydrogen cannot be established here, due to the absence of a well and uniquely defined 
classical Hamiltonian. In particular, the atomic core destroys the symmetry characteristic for 
the hydrogen atom and the Runge-Lenz vector is no more a constant of motion. 

Indeed, experimental data systematically suggest strongly enhanced ionization rates of 
nonhydrogenic (i.e., low angular momentum) alkali Rydberg states as compared to atomic 
hydrogen ||, H, 0, ^, 0, though they also exhibit qualitatively similar features, e.g. of the 
dependence of the ionization yield on the principal quantum number of the atomic state the 
atoms are initially prepared in 10 . On the other hand, a direct comparison of available 



hydrogen and alkali data is somewhat questionable, since relevant experimental parameters 
such as the interaction time of the atom with the field are typically different for different 
experiments. Furthermore, a rigourous theoretical treatment of alkali atoms exposed to 
microwave fields was not accomplished until now. 

It is the purpose of the present letter to outline such a rigourous treatment which allows for 
the first time for a direct comparison of hydrogen and alkali ionization dynamics under precisely 
the same conditions, without adjustable parameters. First results of our numerical experiments 
directly address the above question of quantum-classical correspondence for periodically driven 
alkali atoms. 

Theory. - Let us start with the nonrelativistic Hamiltonian of a one-electron atom exposed 
to a linearly polarized microwave field of (constant) amplitude F and frequency uj, in length 
gauge, employing the dipole approximation and atomic units: 

H{t) = y + Kitom(r) + Fzcosiut, r > 0. (1) 

As this Hamiltonian is periodic in time, we can use the Floquet theorem ||ll|] to find the 
eigenstates ("dressed states") of the atom in the field. After integration over the sohd angle 
we have to solve the time- independent, radial eigenvalue equation 

+ ^^7^ + 2Ktom (r) - 2kuj - 2e^ 1*^;,) 
+FrA,^, + ,)) + FrA, + \^^^+/_^)) = 0, 

with Ae^J M^-i ' ^ = 0' 1' 2, . . . ; fc = ~oo, +oo. (2) 

The additional quantum number k counts the number of photons that are exchanged between 
the atom and the field, and e denotes the quasi-energy of the dressed state 

= 5Zexp(-ifcwt)|*^) = ^exp(-ifcwt)rf,„,(0,0)|^',^^,)/r, (3) 

k k,e 

with Yi „i{9, (p) the spherical harmonics, m denotes the angular momentum projection on the 
field polarization axis and remains a good quantum number, due to the rotational symmetry of 
our problem around the field axis. For all numerical results presented hereafter, its value was 
fixed to m = 0. As immediately obvious from the nondiagonal part of eq. (|^), the interaction 
with the linearly polarised microwave field conserves the generalised parity H = (—1)*'"''^. This 
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just expresses the angular momentum transfer associated with the absorption (emission) of a 
photon. 

As a unique one-particle potential Vatom('') for alkali atoms is unknown, we use a variant 
[|2| of R-matrix theory to describe the interaction of the outer electron with the atomic core. 
Configuration space is divided in two regions: In the internal region, < r < a, the external 
field is negligible compared to the field created by the atomic core, and the details of the 
interaction are unknown. With the help of quantum defect theory the solution of eq. @) 
at r = a can be written as a linear combination of regular and irregular Coulomb-functions 
se,Eir) and Q,_B(r), 

Fe^E{r) = cos{'KSe)se,Eir) + s'm{n S e) eg ^E{r), r = a, (4) 

where the Se are the quantum defects known from spectroscopic experimental data [Q. 
In the outer region, r > a, the difference between the actual atomic potential T4tom(f) and the 
Coulomb potential — 1/r can be neglected. However, the operator <P /dr^ is no more hermitian 
in the reduced range a < r < oo. To overcome this problem, a surface term d{r — a){-§p + Ci) 
is added jlj, |l^ to the diagonal part of (^) . The matching condition between inner and outer 
region at r = a is incorporated in the constant C'l by defining 

C, = (F,,,+fc^(r))(-i)|-F,,,+fe^(r). (5) 
or 

Note that the function Fi E{f) in eq. (^) has to be evaluated at the energy e + ku) in (H), i.e. 
at different energies for different photon indices k. This generalizes the approach outlined in 
[ p^ to periodically driven systems. 

Finally, due to the continuum coupling induced by the external field, all atomic bound states 
turn into resonances with finite ionization rates Fg. In order to extract the latter together with 
the energies e of the atom in the field, we use the method of complex scaling 0. After 
this nonunitary transformation the Floquet Hamiltonian amended by the core induced surface 
term is represented by a complex symmetric matrix, with complex eigenvalues e — 
These are obtained by diagonalization of the complex eigenvalue problem in a real Sturmian 
basis, using an efficient implementation of the Lanczos algorithm. Together with the associated 
eigenvectors they provide a complete description of our problem p5[ |. 

Results. - The described theoretical/mmrerical apparatus is now applied to alkali atoms in 
a microwave field. Since we want to identify the core induced effects in the alkali problem as 
compared to the hydrogen spectrum, we use parameter values which have been employed in 
earlier work on microwave driven Rydberg states ^ of hydrogen. To keep the comparison 
as transparent as possible, we focus on a microwave frequency oj — 1.07171794 x 10^* a.u. 
which is nonresonant with the hydrogen level spacing in the vicinity of the atomic initial state 
with principal quantum number no — 23. The field amplitude is fixed to F = 1.072 x 10~^ a.u., 
slightly below the onset of appreciable (chaos-induced ^) ionization of atomic hydrogen Q. 
This choice of parameters defines a near-integrable phase space structure for the classical 
dynamics of driven hydrogen, with an unambiguous signature in the associated quantum 
energies emerging from the no — 23 manifold. The black dots in fig. ^ illustrate the situation: 
The driving field lifts the angular momentum degeneracy of the substates of the manifold, 
which reorganize according to their localization properties in classical phase space Q. Those 
states with maximum angular momentum and spherical symmetry experience the strongest 
field induced ( "ac-" ) shift in energy, whereas those with maximum radial component of the 
Runge-Lenz vector and "A-symmetry" ^ |l^ remain essentially unaffected by the external 
perturbation. Since the low angular momentum states are strongly mixed by the field (to 
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build states with A-symnictry ^ ) , a new (semiclassical) quantum number p 0| is used to 
label the uq substates of the manifold in the field, p is an integer ranging from to tiq — 1, 
and simply counts the number of quanta enclosed by a semiclassical contour integral along 
the equipotential curves of the adiabatic Hamiltonian which generates the slow evolution 
of angular momentum and Runge-Lenz vector of the classical Kepler ellipse under external 
driving Q. The associated eigenstates exhibit spherical symmetry for p = 0...9, and A- 
symmetry for p = 10 . . . 22, respectively [||. Note that low and high p- values correspond to 
negligible ionization rates of the atom in the field, due to the classical stability of the associated 
trajectories under external driving [Q. Actually, the A-states with large p, which quantize a 
classical straight line orbit perpendicular to the field polarization axis, with maximum modulus 
of the Runge-Lenz vector, display the smallest ionization rates Q]. 

In the presence of a non-hydrogenic core, the Runge-Lenz vector is no more a conserved 
quantity and the A-symmetry defining associated eigenstates of the field free atom is 
destroyed. Therefore, no symmetry argument is available to predict a similar (semiclassical) 
organization of the alkali energy levels under external driving, alike the one observed for atomic 
hydrogen |@. 

Nonwithstanding, our results for lithium Rydberg states exposed to precisely the same 
external perturbation as for the hydrogen results clearly show that the symmetry properties 
of the driven Coulomb problem prevail even in the presence of the core. As evident from the 
open triangles in fig. (a), the hydrogenic part of the lithium manifold exhibits globally the 
same (semiclassical) structure as the hydrogen levels. For low values of p (~ ... 9) this is not 
surprising as the associated classical trajectories (large angular momenta) do not probe the 
atomic core 0. However, for large p- values (~ 10 . . . 20), the classical solution of the Coulomb 
problem does impinge on the nucleus and will certainly suffer scattering off the nonhydrogenic 
core. Yet, in the presence of the field, this scattering obviously mixes states of A type only and 
does not affect the overall separation of the spectrum in spherical and A states, as a remnant of 
the classical phase space structure of the driven Coulomb dynamics. Neither does the presence 
of the core appreciably affect the ionization rates of the dressed states, as obvious from fig. |^ 
(b). Only at p = 10 is there a local enhancement of the width (by approx. one order of 
magnitude), due to the near resonant coupling of the state to the nonhydrogenic eigenstate 
originating from |ri = 41, ^ = 0), via a six-photon transition (similarly, a very weak mutliphoton 
coupling slightly enhances the width of the p = 12 state). In the near integrable regime of the 
classical Coulomb dynamics we are considering here it is precisely this kind of multiphoton 
resonances between nonhydrogenic (low £, such that Si ^ 0) states and hydrogenic manifolds 
which provides a channel for enhanced ionization as compared to atomic hydrogen. Note that 
without such a near resonant coupling, the non-hydrogenic states of a given manifold tend 
to be more stable than the hydrogenic ones, as they are highly isolated in the spectrum. As 
an example, for the same field parameters, the lithium no = 23 £ = ((5^=o = 0.399468) and 
I = \ (Si=\ = 0.047263) ||l^ states exhibit ionization rates ~ 10^^^ a.u. as small as the 
most stable substates of the hydrogenic manifold of fig. |l|. A detailed analysis of enhanced 
ionization via core-induced multiphoton resonances will be provided elsewhere. 

Closer inspection of fig. |l| (a) shows additional structure in the alkali spectrum, on top of 
the globally hydrogen-like structure: for large values oip (> 11), the alkali levels are shifted 
with respect to the hydrogenic energies. These shifts can be recovered by diagonalization of 
the hydrogen problem within the restricted subspace spanned by the hydrogenic levels of the 
alkali Rydberg manifold ||l^, |2^, |2l|, In other words, the shifted energies are the solutions 
of the eigenvalue equation 



PHhydP\^''«} = {E + fcow)|$*=°). 



(6) 
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where i?hyd is obtained from from (1) setting Kitom(?') = ^1/^, f g]0, cxd[, and P the projector 
onto the hydrogenic subspace of the alkah manifold labeled by the principal quantum number 
no and the photon number fco. Such a procedure is legitimate as long as the states emerging 
from the nonhydrogenic part of the alkali manifold have vanishing overlap with the complete 
hydrogen manifold emanating from (no,fco). This condition is fulfilled for the driving field 
strength considered here. 

Solving (0) for E is tantamount to finding the roots of 



det(Q 



-ffhyd -{E + koLj) 



(7) 



with Q = 1 — P the projector onto the orthogonal complement of the hydrogenic subspace for 
given {no, fco)- Without loss of generality we choose fco = hereafter. Consequently, for one 
single non- vanishing quantum defect Sig , (Q) becomes 



E 



e-E 



= 0, 



(8) 



where |no,£o) spans the orthogonal complement of the hydrogenic subspace of the alkali atom 
within the (no, fco = 0) manifold. Note that or (^) have to be evaluated for different values 
of the generalized parity 11, and that we have to solve separately for ^o = and = \/vci 
order to recover the level shifts observed for lithium in fig. |^ (the ^o = 2 and i'o = 3 states of 
lithium remain within the range of P, due to their negligible quantum defects 5^=2 — 0.002129 
and 5i=2 = —0.000077 Q, at the given field strength). Fig. ^ (a) shows the result of the 
projection method, compared to the exact numerical result - the agreement is very good. 
Since the low p states essentially exhibit spherical symmetry with large angular momentum 
projection, their overlap with |no,^o = 0(^o = 1)) vanishes and their energies remain unshifted 
as compared to the hydrogen results. 

The scenario which we described for lithium also applies for the heavier alkali elements, as 
illustrated in figs. ^ (b) and (c). Here we plot the shifts of the exact energies of sodium and 
rubidium with respect to the hydrogen levels, as they emerge from the no — 23 manifold, for 
precisely the same field parameters as used for the lithium results. Since for these elements 
also the ti =2 (sodium) and the = 3 (rubidium) states are separated from the hydrogenic 
manifold due to their large quantum defects, the range of Q in (^) is two-dimensional and the 
evaluation of the determinant yields the expression 



E 



^no,£o|*^=°)P V- l(no,^i|*^=°)l' 



E 



■E 



E 



E 



(no,4|^^=°)(^^=°|no,^i) 

e-E 



= 0. (9) 



Again, the solution of (^) gives very good agreement with the numerical result. In addition, 
we note that the larger the dimension of the range of Q, the smaller the values of p for which 
the alkali levels are shifted as compared to the hydrogen energies. This is a consequence of 
the dominance of small t components in large p states and of large € components in small p 
states, since the heavier the element the larger the i values affected by non-negligible quantum 
defects. 



Summary. - In conclusion, the energy levels of alkali Rydberg states emerging from the hy- 
drogenic no-manifold clearly reflect the phase space structure of the microwave driven Coulomb 
problem, despite the presence of a symmetry breaking atomic core. Also the ionization rates 
of the atoms reflect the underlying classical phase space structure, with the exception of local 
enhancements due to multiphoton resonances with nonhydrogenic sublevels of other manifolds. 
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We have checked that the observed structure is robust under changes of the driving field 
amphtude, up to values where adjacent n-manifolds start to overlap. 

We thank Dominique Delande and Ken Taylor for fruitful discussions and an introduction 
to the R-matrix approach of . 
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Fig. 1. - Energies (a) and ionisation rates (b) of Rydberg states of lithium (triangles) and of atomic 
hydrogen (dots) exposed to a linearly polarized microwave field of frequency uj — 1.07171794 x 10^* a. u. 
and amplitude F — 1.072 x lO^'^a.u., for principal quantum number no — 23 and angular momentum 
projection m = on the field polarization axis. The lithium spectrum lacks two of the 23 substates 
of the manifold, due to the quantum defects 5t=o = 0.399468 and 5i^i — 0.047263 of the 1 = 
and £ — 1 states, respectively. The quantum defects 6e-2 = 0.002129 and 5^=3 = —0.000077 
are negligible compared to the field induced splitting of the no = 23 manifold (field-free energy 
E23 — —9.452 X 10"* a.u.). Both specti^a almost coincide (in energy and ionisation rate) even for larger 
values (p > 10) of the (semiclassical M) quantum number p, despite the fact that the localization 
properties of the associated eigenstates (close to the plane defined bv the field polarization axis) 
originate in the dynamical symmetry of the — 1/r Coulomb potential [jl^]. The latter is destroyed by 
the presence of a nonhydrogenic core in alkali atoms. The ionization rate of the p = 10 state of lithium 
is locally enhanced by approx. one order of magnitude with respect to the corresponding hydrogen 
eigenstate, due to a six-photon resonance with the |n = 41, ^ = 0) state. 



Fig. 2. - Shifts Saik — ^^hyd of the energies i^aik of lithium (a, triangles), sodium (b, diamonds), and 
rubidium (c, squares) as compared to those, -Ehyd, of the no = 23 manifold of atomic hydrogen in a 
linearly polarized microwave field, with the same parameters as in fig. |l|. Quantum defects employed 
for the sodium results: 5e=o = 1.347964, Se=i = 0.85538, Si=2 = 0.015543, = 0.001453, and for 
rubidium; 5«=o = 3.1311, 5^=1 = 2.6415, 5^=2 = 1.3472, de=3 = 0.016312 ||lj. Consequently, three 
respectively four energy levels are missing in (b) and (c). The nonvanishing shifts for large p > 9 values 
can be accounted for by projecting out the low £ components (i.e. the ones with core induced energy 
shifts large with respect to the field induced splitting of the no — 23 manifold) of the no-manifold, as 
indicated by the crosses, see eqs. (8) and (9) . The agreement between this perturbative approach and 
the exact quantum results is always better than the average level spacing of the hydrogen manifold 
(dots in fig. 1), except for the relatively large discrepancy at p = 11, in (c). The latter is due to a 
multiphoton resonance between the alkali eigenstate and a nonhydrogenic (low £) state. 



